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PREFACE 


SARD,  INC.  the  research  and  development  subsidiary  of  6ATX,  has  prepared 
this  report  for  the  Federal  Emergency  Management  Agency  (FEMA)  in  Washington, 
D.C.  Originally,  the  study  was  contracted  for  by  the  Defense  Civil  Prepared¬ 
ness  Agency  (DCPA),  which  has  since  become  part  of  FEMA.  Mr.  Donald  Bettge 
of  FEMA  served  as  Project  Monitor  during  the  entire  program. 

This  report  describes  the  results  of  the  experimental  and  analytical 
studies  that  were  conducted  to  aid  in  the  preparation  of  a  methodology  to 
determine  the  potential  of  natural  ventilation  to  ventilate  upgraded  shelters. 
Scale  model  experimental  tests  were  conducted  using  a  unique  low  speed  wind 
tunnel  that  was  designed  and  built  especially  for  this  program.  Experimental 
data  was  then  incorporated  into  a  mathematical  model  which  can  be  used  to  assess 
the  adequacy  of  natural  ventilation  to  limit  the  internal  environment  of 
certain  simple  upgraded  shelter  configurations  to  acceptable  limits. 

Individuals  at  CARD  which  participated  in  this  program  include: 

R.  H.  Henninger  -  Project  Engineer 

Dr.  R.  J.  Tsai  -  Technology  Review  and  Mathematical  Modeling 

Dr.  S.  F.  Fields  -  Experimental  Modeling 

Dr.  C.  K.  Krishnakumar  -  Experimental  Modeling 

J.  B.  Koh  -  Data  Reduction 

R.  J.  Honegger  -  Wind  Tunnel  Design 

GARD  wishes  to  thank  Mr.  Bettge  and  FEMA  for  the  opportunity  to  have 
undertaken  this  study. 


Director,  Contract  Programs 


Respectfully  submitted, 

R.  H.  Henninger 
Project  Engineer 


ABSTRACT 


Experimental  and  analytical  investigations  were  conducted  during  this  study 
for  the  purpose  of  evaluating  the  adequacy  of  natural  ventilation  in  upgraded 
shelters.  A  unique  low-speed  wind  tunnel  which  uses  photographic  measurement 
techniques  for  flow  tracing  of  neutrally  buoyant  bubbles  through  openings  was 
utilized  to  conduct  scaled  model  tests  of  three  shelter  models  to  determine 
the  ventilation  air  throughput  CFM  as  a  function  of  wind  speed,  relative  wind 
approach  angle  and  opening  pattern.  The  results,  as  summarized  below,  indicate 
that  sizeable  ventilation  rates  are  achieveable  at  low  wind  speeds  and  based 
upon  data  reported  by  FEMA  would  result  in  adequate  ventilation  rates  to 
meet  the  82°F  effective  temperature  and  90%  adequacy  criteria  for  all  but  the 
southeast  portion  of  the  U.S.  Further  research  is  required  however,  to  determine 
the  effects  of  internal  flow  resistances  provided  by  partition  walls  and 
occupants,  air  stratification,  areas  and  locations  of  openings,  etc.  before 
these  results  can  be  accepted  with  confidence. 


ESTIMATED  VENTILATION  FLOWS  RATES  ACHIEVABLE  WITH  NATURAL 
VENTILATION  BASED  ON  SCALED  MODEL  TESTS  FOR  WIND  SPEEDS  IN 
EXCESS  OF  5  mph  AMD  OCCUPANT  DENSITY  OF  10  SQ.  FT.  PER  PERSON 
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FEMA  indicates  that  7.5  to  40  CFM  per  person  is  adequate  to  maintain  82°F 
effective  temperature  and  90%  adequacy  (Ref:  ASHRAE  Applications  Handbook, 
1978,  Chapter  12,  Figure  13)  GAHO,  INC. 
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Section  1 


INTRODUCTION 


1.1  Background 

The  survival  of  the  occupants  of  a  shelter  is  dependent  upon  the 
following  environmental  factors: 

1)  the  nuclear  radiation  level  in  the  shelter, 

2)  the  chemical  composition  of  the  air  in  the  shelter,  and 

3)  the  thermal  level  in  the  shelter. 

Among  these  three  requirements,  the  second  and  the  third  are  related  to  the 
ventilation  rate  of  the  shelter.  Previous  studies  by  the  Defense  Civil 
Preparedness  Agency  (now  part  of  the  Federal  Emergency  Management  Agency) 
indicated  that  a  ventilation  rate  of  approximately  3  CFM  per  occupant 
would  suffice  for  chemical  balance.  Since  this  ventilation  rate  is  easily 
attainable,  the  chemical  composition  of  the  shelter  air  is  adequately  pro¬ 
vided  for.  Thus,  the  main  environmental  control  problem  in  a  shelter  is 
the  provision  of  sufficient  ventilation  to  protect  the  occupants  against 
excessive  thermal  levels  within  the  shelter. 

Ventilating  a  shelter  with  ambient  air  to  control  the  thermal  level 
within  the  shelter  might  be  accomplished  in  two  ways: 

1)  Natural  Ventilation  -  ventilation  through  open  windows,  doors 
stairwells,  etc.  caused  by  atmospheric  winds  or  thermal  effects 

2)  Forced  Ventilation  -  ventilation  by  artifical  means  using 
mechanical  devices. 

If  natural  ventilation  can  be  shown  to  be  reliable,  one  approach  might 
be  to  utilize  Pedal  Ventilator  Kits  (PVK)  and  Kearny  Pump  Kits  (KPK) 

(Ref.  1.1)  only  in  those  locations  and  shelters  where  natural  ventilation 
would  not  be  adequate.  The  need  arose  therefore,  for  an  accurate  assess¬ 
ment  technique  to  be  developed  which  could  determine  the  adequacy  or 
non-adequacy  of  natural  ventilation  for  upgraded  shelters,  shelters 
whose  ventilation  characteristics  had  not  been  previously  studied. 
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1.2  Study  Objective  and  Scope 

Previous  studies  concerning  natural  ventilation  effects  of  fallout 
shelters  have  dealt  with  typical  above  and/or  below-ground  fallout  shelters 
in  existing  structures;  this  study  deals  with  upgraded  shelters,  i.e.  shelters 
having  full  earth  berms  with  at  least  6  inches  of  earth  overhead.  It  is 
questionable  if  previous  field  test  data  (Section  2)  obtained  for  sharp- 
edged  structures  are  applicable  to  the  sloped-sided  configuration  of  upgraded 
shelters.  The  objective  of  this  study  was  therefore  to  experimentally  and 
analytically  assess  the  potential  of  natural  ventilation  in  upgraded  shelters. 

1.3  Study  Approach 

The  wind  induced  ventilation  air  flow  rate  through  a  building  is  a 
complex  function  of  several  geometric  and  flow  variables.  Important  among 
these  variables  are  wind  speed,  wind  direction,  boundary  layer  velocity 
profile  of  approaching  wind,  building  geometry,  areas  and  locations  of  windows 
and  doors,  internal  obstacles  within  the  building  and  the  nature  and  proximity 
of  neighboring  buildings  and  obstructions.  Where  the  building  is  situated  in 
open  country,  the  problem  simplifies  to  some  extent.  In  that  case,  reasonable 
assumptions  can  be  made  for  the  boundary  layer  profile  of  the  approaching  wind. 
Further,  if  the  building  is  sharp  edged,  the  flow  separation  points  are  well 
defined.  This  is  not  true  of  flows  over  buildings  with  curved  exteriors  like 
the  bermed  buildings  of  the  present  study.  Even  for  buildings  of  relatively 
simple  geometries,  however,  the  velocity  and  pressure  profiles  in  and  around 
these  buildings  are  so  complex  as  to  render  a  complete  analytical  or  numerical 
solution  impractical. 

A  common  approach  used  to  study  problems  of  this  nature  is  to  construct 
a  theoretical  model  and  complement  it  with  data  from  carefully  conducted  exper¬ 
iments.  Data  acquisition  from  experiments  of  full  size  buildings  alone  is 
prohibitively  expensive  and  time  consuming.  The  solution  therefore  lies  in 
conducting  tests  on  properly  designed  scaled  model  systems,  feeding  the  data 
into  a  theoretical  model  and  finally  testing  the  results  with  representative 
field  data.  Often,  feed  back  between  model  and  field  tests  may  be  necessary 
to  correct  inaccuracies  in  the  model  system  or  to  interpret  field  data.  Such 
an  approach  has  been  adapted  here  where  a  series  of  tests  were  designed  to 
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yield  a  clear  understanding  of  the  effects  of  the  important  variables  individ¬ 
ually,  and  in  groups,  on  ventilation  air  flow  rates.  The  building  geometry  was 
kept  relatively  simple  and  internal  flow  resistances  were  excluded.  Natural 
ventilation  due  to  thermal  effects  are  also  being  ignored  for  the  present. 
Information  gained  from  these  studies  together  with  data  from  future  tests 
on  building  models  with  internal  resistances  should  enable  reasonably  accurate 
predictions  to  be  made  of  ventilation  air  flow  rates  through  full  size 
buildings  of  the  types  considered. 


1.4  Study  Plan 

The  study  plan  employed  for  this  project  was  structured  around  the  exper¬ 
imental  and  analytical  approach  described  in  the  previous  section  and  included 
the  following  work  tasks: 


Task  1  - 
Task  2  - 

Task  3  - 


Task  4 


Task  5 
Task  6 


Technology  review 

Identify  applicable  analytical  flow  models,  flow  parameters 

and  shelter  configurations 

Scaled  model  testing 

3A  -  Formulate  test  plan 

3B  -  Investigate  scaling  effects 

3C  -  Design  and  build  test  stand  and  models 

3D  -  Calibrate  wind  tunnel 

3E  -  Conduct  tests 

3F  -  Reduce  data  and  interpret  results 
Analytical  model  testing 

4A  -  Develop  analytical  natural  ventilation  assessment  model 
4B  -  Incorporate  experimental  data  into  model 
Assess  adequacy  of  natural  ventilation  based  on  results 
Prepare  final  report  and  recommendations. 
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Section  2 


TECHNOLOGY  REVIEW 


Natural  ventilation  as  applied  to  buildings  (Refs.  2.6  through  2.14) 
found  widespread  application  back  in  the  1920-40  era  where  ventilation  due 
to  wind  and  thermal  effects  were  routinely  incorporated  into  building  designs. 

The  acceptance  and  general  use  of  air  conditioning  since  that  period 
however,  has  limited  its  use  to  only  industrial  buildings.  Presently,  there 
seems  to  be  a  resurgence  of  interest  in  natural  ventilation  as  a  building 
energy  conservation  technique. 

The  state  of  the  art  of  natural  ventilation  calculation  techniques  is 
described  in  the  1977  ASHRAE  Handbook  of  Fundamentals.  These  procedures  have 
their  roots  in  data  and  experience  gained  from  the  1920 ' s  and  30 's  and  have 
changed  little  since  that  time.  Research  in  the  last  decade  has  concentrated 
in  an  area  directly  related  to  natural  ventilation,  that  being  building  infil¬ 
tration.  Both  are  wind  induced  and  both  deal  with  ambient  air  movement  through 
aperatures.  Natural  ventilation  is  generally  controllable  and  desireable, 
whereas  infiltration  cannot  be  effectively  controlled  and  is  undesireable.  A 
survey  was  made  of  infiltration  studies  to  determine  the  applicability  of  results 
and  mathematical  models  to  this  present  study.  Mathematical  models  developed 
as  part  of  these  studies  are  varied  and  generally  can  be  applied  to  only  certain 
building  types.  This  section  summarizes  those  models  which  have  been  developed 
for  treating  both  natural  ventilation  and  infiltration  in  buildings  and  structures. 

1 

I 

2.1  ASHRAE  Natural  Ventilation  Model  And  Application  to  Above-Ground  Fallout  ! 

Shelters  i 

it 

Extensive  experimental  and  analytical  studies  of  natural  ventilation  for 
above-ground  fallout  shelters  were  conducted  by  the  Defense  Civil  Preparedness 
Agency  back  in  the  1960's.  These  studies  utilized  a  relationship  obtained  from 
the  continuity  equation  to  determine  wind-induced  air  flow  through  ventilation 
openings.  This  relation,  based  upon  the  presentation  in  the  1965  ASHRAE 
Guide  and  Data  Book,  is  given  below  and  still  represents  the  state  of  the  art 
as  is  witnessed  in  the  1977  ASHRAE  Handbook  of  Fundamentals  (Ref.  2.1). 
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Q 

where  Q 

B 


(1  +  B)  EAV 

volume  flow  rate  (ft  /mln/occupant) 

fractional  multiplying  factor  which  arises  from  unequal  Inlet 
and  outlet  areas 

effectiveness  factor 

Atot^N  2 

smaller  of  the  air  Inlet  area  or  air  outlet  area  (ft*) 

number  of  occupants 


V  *  atmospheric  wind  speed  ( ft/ml n). 


The  value  of  B  is  dependent  upon  the  ratio  of  inlet  area  to  outlet  area,  or 
vice  versa,  whichever  is  greater  that  one.  This  dependence  Is  not  linear,  but 
is  as  Indicated  in  Figure  2.1.  Increasing  inlet  area  over  outlet  area,  or 
vice  versa,  will  increase  the  Induced  ventilation,  but  not  in  proportion  to 
the  added  area. 


Figure  2.1  INCREASE  IN  FLOW  CAUSED  BY  INEQUALITY  OF  INLET  AND  OUTLET  AREAS 

The  E- factor  in  the  above  equation  represents  the  effectiveness  of 
openings  and  ASHRAE  Indicates  that  E  should  be  taken  as  0.50  to  0.60  for  winds 
perpendicular  to  the  openings,  and  0.25  to  0.35  for  diagonal  winds.  Experimental 
tests  conducted  by  DCPA  for  4  different  shelters  (Refs.  2.2  through  2.5) 

Indicated  that  as  atmospheric  wind  speed  Increased,  the  air  flow  induced  by  the 
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atmospheric  wind  did  not  follow  the  expected  linear  relationship.  This  can  be 
seen,  for  example,  by  examining  Figure  2.2  which  presents  the  results  for  the 
Evanston,  Illinois  test  (Ref  2.4).  In  this  rectangularly-shaped  aboveground 
shelter,  all  windows  had  the  same  openable  area.  In  Test  1,  four  windows  were 
open  (one  on  each  wall)  and  in  Tests  2  and  3,  two  windows  were  open  (one  on  each 
opposite  wall  with  different  opposite  walls  in  Tests  2  and  3).  The  test 
results  show  a  scattered  distribution  and  departure  from  linearity,  especially 
at  higher  atmospheric  wind  speeds.  This  behavior  Is  probably  <*je  to  several 
other  effects  that  are  present  but  not  accounted  for  In  the  ASHRAE  equation, 
e.g. , 

1)  the  effect  of  the  atmospheric  wind  direction, 

2)  the  type  of  window  openings, 

3)  the  Interior  structural  elements, 

4)  the  exterior  structures  that  disturb  the  surface  wind  velocity 
field,  and 

5)  the  variation  In  height  between  the  shelter  location  and  the 
location  at  which  the  atmospheric  wind  speed  and  direction  are 
determl ned. 

2.2  Other  Natural  Ventilation  And  Building  Infiltration  Flow  Models 

Some  100  or  more  technical  papers  relating  to  natural  ventilation  and 
Infiltration  of  buildings  and  structures  were  reviewed  to  determine  the  scope 
of  knowledge  related  to  these  subjects.  A  compilation  of  the  more  pertinent 
papers  are  presented  in  the  reference  list  found  at  the  end  of  this  section. 
Appendix  A  sunmarizes  the  mathematical  flow  models  that  were  gleaned  from  this 
review. 

Since  the  1960's,  little  work  has  been  carried  on  In  the  specific  area  of 
natural  ventilation  of  buildings  except  for  the  four  DCPA  studies  previously 
referenced.  Recent  studies  have  concentrated  rather  In  the  related  area  of 
infiltration  of  ambient  air  Into  residential  and  commercial  buildings  through 
their  exterior  envelopes.  Ross  and  Grimsrud  (Ref.  2.53)  compiled  a  thorough 
sunmary  of  past  and  recent  works  conducted  In  this  area  and  should  be  consulted 
for  a  more  complete  discussion.  Many  of  the  Infiltration  models  developed  as 
part  of  those  studies  were  specific  to  a  particular  building  or  building 
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classification  and  as  such  are  not  applicable  to  our  study. 

The  ASHRAE  linear  model  discussed  previously  represents  the  simplest  of 
the  flow  models  In  existence.  It  assumes  that  the  air  flow  through  a  building 
Is  a  linear  function  of  wind  speed  only,  leaving  all  other  effects,  such  as 
wind  direction,  type  of  opening,  ratio  of  opening  area  to  wall  area,  etc.,  to 
be  accounted  for  by  the  Effectiveness  Factor  (E). 

In  1955,  G.  A.  Mackslmov  (Ref.  2.13)  formulated  an  equation  for  air  flow 
through  buildings  as  a  function  of  static  pressure  difference  between  the 
windward  and  leeward  sides  of  the  building.  This  was  developed  through 
application  of  Bernoulli's  Law  and  resulted  in  a  non-linear  function  of  air 
flow  to  wind  speed. 

Many  authors  have  attempted  to  develop  specific  non-linear  equations  by 
using  results  of  experimental  field  tests  (Refs.  2.7,  2.13,  2.19,  2.20,  2.38, 
2.46,  2.47,  2.50,  2.51,  2.54,  2.56,  2.57).  Usually  the  dynamic  flow 
coefficients  obtained  were  for  specific  building  configurations  and  character¬ 
istics  and  could  not  be  generalized  to  other  types  of  buildings  and  conditions. 

Although  most  authors  recognize  that  air  flow  through  a  building  Is  a 
function  of  wind  direction,  few  have  Included  this  effect  directly  as  an  inde¬ 
pendent  variable.  Some  (Refs.  2.6,  2.21,  2.25,  2.46,  2.50,  2.56)  have 
included  the  relative  angle  between  the  wind  direction  and  the  normal  to  the 
wall  in  linear  or  non-linear  equations  but  have  developed  pressure  coefficients 
for  only  one  reference  direction.  Pressure  coefficients  are  of  course  a 
function  also  of  building  configuration,  placement  of  openings,  ratio  of  wall 
and  opening  areas,  etc.  These  have  usually  been  neglected. 

More  general  equations  for  residential  application  have  been  developed 
recently  to  account  for  such  other  effects  as  frequency  of  doorlng  opening, 
presence  of  fireplaces  and  chimneys,  operation  of  ventilators  and  exhaust  fans, 
furnace  operation,  etc.  (Refs  2.10,  2.19,  2.20,  2.21,  2.28,  2.38,  2.46,  2.48, 
2.50,  2.51).  A  good  example  of  this  approach  was  that  done  by  Ohio  State 
University  (Ref.  2.46)  for  residential  application.  As  Is  true  for  this  work 
as  well  as  the  others,  the  coefficients  determined  by  field  tests  for  Insertion 
Into  the  polynomial  expressions  are  specific  to  the  structures  studied  and 
especially  to  the  quality  of  construction. 


Computerized  models  developed  by  Konstantinova  (Ref.  2.20)  and  Sander  and 
Tamura  (Ref.  2.37  and  2.38)  using  numerical  methods  allowed  for  a  more  detailed 
solution  to  the  problem  by  being  able  to  handle  more  complex  structures  where 
a  series  of  equations  must  be  solved  simultaneously.  But  still  these  models 
required  that  pressure  coefficients  be  determined  independently  as  a  function 
of  wind  speed  and  direction.  Computerized  models  also  gave  the  opportunity 
to  utilize  recorded  hourly  weather  data  and  to  study  the  Independent  effects 
of  certain  variables  as  they  are  varied  over  their  range  of  occurrence. 

As  will  be  seen  In  Section  3,  the  type  of  structure  under  study  In  this 
program,  i.e.,  bermed  one-story  buildings,  are  not  like  the  sharp-edged 
structures  studied  to  date  but  rather  have  curved  or  sloped  sides  which  give 
them  flow  characteristics  different  from  typical  buildings.  Dynamic  pressure 
coefficients  and  coefficients  determined  for  polynomial  expressions  previously 
discussed  therefore  made  application  of  these  data  to  our  problem  questionable. 
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Section  3 


SCALED  MODEL  TESTING  PROGRAM  AND  RESULTS 


3.1  Approach 

The  model  tests  in  this  study  were  designed  with  the  object  of  gaining 
insight  into  the  phenomena  of  wind  flow  through  and  around  bermed  shelter 
buildings  and  formulating  relationships  between  ventilation  flow  rate  and 
major  parameters  affecting  it.  It  was  realized  at  the  outset  that  measure¬ 
ment  of  volume  flow  rates  through  the  model  openings  would  be  difficult. 

Flow  patterns  in  the  vicinity  of  the  openings  are  very  complex,  flow  areas 
are  not  well  defined,  and  conventional  probes  are  not  desirable.  In  view  of 
these  problems,  the  goal  for  this  initial  study  was  to  develop  a  viable 
measurement  technique  and  apply  it  to  simple  models.  This  has  been  success¬ 
fully  achieved  in  this  study. 

The  technique  developed  by  GARD  consists  of  tracing  the  flow  of  neutrally 
buoyant  bubbles  through  models  and  recording  their  path  lines  using  motion 
photography.  These  recordings  can  be  analyzed  using  a  variable  speed 
projector  to  yield  valuable  quantitative  and  qualitative  data. 

3.2  Scaling  Considerations 

The  basic  requirement  in  model  testing  is  dynamic  similarity  between  the 
full-scale  and  model  systems.  Dynamic  similarity,which  implies  similarity  of 
geometries  and  flow  fields, is  a  prerequisite  for  translating  model  data  to 
the  full-scale.  In  most  cases,  perfect  similarity  is  impractical.  However, 
it  often  suffices  to  ensure  similarities  only  with  regard  to  the  dominant 
characteristics  of  the  two  flow  systems.  A  dimensional  analysis  of  the 
problem  shows  that  dynamic  similarity  can  be  achieved  by  ensuring  that  1) 
the  model  and  full-scale  systems  are  geometrical ly  similar,  2)  values  of 
Reynolds  numbers  based  on  building  height,  width  and  length  are  the  same  for 
the  model  and  full-scale  systems,  3)  velocity  profiles  of  the  approach  wind 
boundary  layers  are  similar  (which  means  that  the  exponent  of  the  power-law 
distribution  is  the  same  for  the  two  systems),  and  4)  ratios  of  building 
height  to  wind  boundary  layer  thickness  are  the  same  in  the  two  systems. 

(In  open  terrain  or  bush  country,  the  wind  boundary  layer  is  about  60  to  89 
feet  thick,  which  is  4  to  5  times  the  height  of  a  typical  single-story 
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building.)  Equality  of  Reynolds  numbers  is  crucial  only  in  the  laminar  and 
transitional  regimes  of  flow.  In  these  regimes,  the  flow  structure  behind 
edges  and  corners  is  very  sensitive  to  Reynolds  number.  Sharp-edged  buildings 
offer  a  simplification  in  that  the  flow  separation  points  are  well  defined. 

This  is  not  true  of  bodies  with  curved  edges  like  the  bermed  shelters  of  the 
present  study.  The  strong  influence  of  Reynolds  number  on  wake  flow  patterns 
and  separation  points  associated  with  curved  surfaces  is  well  illustrated  in 
the  experimental  plots  of  drag  coefficient  versus  Reynolds  number  for  cylinders 
(Ref.  3.1  and  3.2).  When  the  approach  stream  is  in  the  fully  turbulent  regime, 
the  drag  coefficient  is  nearly  independent  of  Reynolds  number,  (Ref  3.2).  This 
is  because  the  major  component  of  drag  -  the  contribution  from  wake  flow  - 
approaches  an  asymptotic  value  at  those  Reynolds  numbers.  Since,  the  earth's 
boundary  layer  is  always  turbulent,  it  suffices  in  the  present  study,  to 
ensure  that  model  Reynolds  numbers  are  large  enough  to  guarantee  good  turbulent 
mixing  in  the  boundary  layers,  thereby  negating  the  possibility  of  regular 
vortex  shedding  of  the  Karman  vortex  street  type  (Karman  vortex  street  behind 
cylinders  are  observed  only  at  Reynolds  numbers  of  5,000  or  less). 

Two  different  model  systems  -  a  wind  tunnel  (or  air  model)  and  a  water 
tunnel  -  were  considered  as  possible  choices.  Each  system  has  some  advantages 
and  disadvantages  with  respect  to  the  other.  For  example,  pressure  and  velocity 
measurements  can  be  made  more  accurately  in  a  water  tunnel.  Flow  visualization 
is  also  very  effective  in  water.  However,  major  considerations  like  tunnel 
blockage,  wall  boundary  layer  interference  effects,  problems  of  handling  large- 
flow-rate  hydraulic  circuits  and  generating  thick  boundary  layers,  favored  the 
selection  of  a  wind  tunnel  over  a  water  tunnel. 

3-3  Description  of  Low  Speed  Wind  Tunnel 

GARD's  wind  tunnel  was  designed  specially  to  suit  air  flow  studies 
involving  simulation  of  wind  velocity  profiles.  The  following  considerations 
strongly  influenced  its  design. 

1)  Tunnel  width  and  height  must  be  large  enough  to  accommodate  large 
models  without  creating  significant  blockage  or  boundary  layer 
interference  effects  (large  models  are  necessary  to  permit 
reasonably  accurate  and  meaningful  measurements).  If  the  blockage 
of  the  flow  area  caused  by  the  model  is  more  that  about  50/-  of  the 
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tunnel's  cross-section,  velocities  significantly  greater  than  the 
approach  wind  free  stream  velocity  will  be  generated  in  the  test 
section  (Ref.  3.3).  Interaction  of  the  stream  lines  deflected 
off  the  edges  of  the  model  with  the  tunnel  wall  boundary  layers 
could  also  lead  to  errors. 

2)  Thick  boundary  layers  simulating  the  earth's  wind  boundary  layer 
should  be  obtainable  in  the  tunnel's  test  section  (typical  values 
for  the  thickness  of  the  earth's  wind  boundary  layer  in  open 
terrain  or  bush  country  range  from  60  feet  to  80  feet). 

3)  The  tunnel  should  be  readily  adaptable  to  visualization  techniques. 

GARD's  low  speed  wind  tunnel  is  situated  in  a  very  large  room  which  forms 
the  return  circuit  between  the  inlet  and  the  exhaust.  A  photograph  and  a 
schematic  of  the  wind  tunnel  are  shown  in  Figures  3.1  and  3.2,  respectively. 
The  tunnel  consists  of  a  convergent  intake  section,  4  rectangular  sections 
(60"  width  x  30"  height  x  96"  length),  a  transition  section  at  the  exhaust 
end,  and  a  33"  diameter  centrifugal  blower  which  operates  in  the  suction  mode 
and  exhausts  the  air  through  a  90°  elbow.  Free  stream  velocities  of  up  to  20 
fps  can  be  generated  in  the  test  section  with  the  present  blower. 

The  convergent  intake  includes  a  honey-comb  entrance  panel  followed  by  a  fine 
mesh  screen.  This  arrangement  straightens  the  flow  and  makes  the  turbulence 
intensity  uniform.  A  short  distance  downstream  of  the  intake  section  is  a 
manifold  that  blows  small  jets  of  air  counter  to  the  main  stream  at  a  downward 
angle.  The  succeeding  sections  have  square  bars  laid  on  the  floor  at  regular 
intervals  spanning  the  entire  width.  The  counter  jets,  together  with  the 
lateral  bars,  have  been  found  to  efficiently  generate  and  sustain  thick 
boundary  layers  of  the  magnitude  required  in  the  present  study  for  wind 
velocity  profile  simulation.  This  system,  also  called  the  Momentum  Defect 
Generating  System,  was  instituted  after  a  careful  analysis  of  several  tech¬ 
niques  that  have  been  tried  in  the  past  for  the  same  purpose  (Ref.  3.4,  3.5 
and  3.6).  The  test  section  is  fabricated  from  clear  Plexiglas  sheets  and  has 
an  access  door  located  on  the  side  wall.  Sections  1  and  4  are  clear  on  the 
sides  and  Section  2  is  clear  on  the  top  and  the  sides.  The  transition  section, 
following  Section  4,  consists  of  a  diffuser,  a  vibration  isolator,  a  large 
area  chamber  with  a  honey-comb  panel  and  an  adapter  leading  to  the  blower. 
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Figure  3.1  LOW  SPEED  WIND  TUNNEL  LOCATED  IN  GARD  LABORATORY  BUILDING,  NILES,  ILLINOIS 


Figure  3  .2  SCHEMATIC  OF  WIND  TUNNEL 


A  damper  fixed  to  the  elbow  on  the  blower  exhaust  regulates  the  tunnel  through¬ 
put  flow  rate.  Figures  3.3  through  3. 10  provide  a  walk-through  of  the  wind 
tunnel  and  its  instrumentation. 

Instrumentation  of  the  tunnel  consists  of  the  following  : 

1)  Pitot-static  tubes  P1  and  P2  (Figures  3.2  and  3-7).  These  probes 
determine  the  pressure  and  velocity  profiles  in  the  tunnel's 
cross-section  upstream  of  the  model  and  in  its  vicinity.  They  are 
connected  to  a  precision  manometer  that  can  accurately  read  pressure 
differentials  to  0.002  inches  of  water.  Probe  Pj  can  traverse  the 
entire  width  and  height  of  the  tunnel's  cross-section  while  probe  P2 
can  additionally  traverse  a  distance  of  4*s  feet  in  the  longitudinal 
direction. 

2)  Bubble  generator  unit  (Figures  3.8,  3.9  and  3.11).  This  unit 
delivers  neutrally  buoyant  bubbles  of  high  stability  through  two 
bubble  release  tubes  located  about  30  inches  upstream  of  the  model. 

The  unit  consists  of  high  pressure  cylinders  of  helium  and  air,  a 
3-channel  bubble  generator  console,  two  bubble  heads,  vortex  chambers 
and  bubble  release  tubes.  Within  each  bubble  head,  helium  passes 
through  a  central  hypodermic  tube.  Helium  also  forces  a  bubble  film 
solution  stored  in  the  console  through  an  annulus  surrounding  the 
hypodermic  tube  to  form  helium-filled  bubbles.  The  bubbles  are  blown 
off  the  tip  by  compressed  air  passing  through  the  outermost  annular 
passage  in  the  head.  By  properly  metering  the  helium,  bubble  film 
solution  and  compressed  air  flow  rates,  and  setting  the  helium  and 
air  supply  pressures,  neutrally  buoyant  bubbles  of  nearly  uniform 
size,  from  1/16  to  1/4  inches  in  diameter,  can  be  continuously 
generated  from  each  head.  Lighter  and  heavier  bubbles  are  filtered 
off  by  the  vortex  chambers,  so  that  only  nearly  neutrally  buoyant 
bubbles  reach  the  release  tubes.  These  bubbles  are  much  more  stable 
than  normal  soap  bubbles.  Although,  some  of  the  bubbles  do  burst  on 
contacting  solid  boundaries,  a  good  many  of  them  successfully  trace 
the  flow  completely  through  and  around  the  model.  The  bubbles  have 
good  reflectivity,  so  that  with  proper  lighting  they  can  be  photo¬ 
graphically  recorded. 
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Figure  3.3  HONEY-COMB  FLOW  STRAIGHTENER 


CONVERGENT  INTAKE  CONE 
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Figure  3.5  COUNTER-JET  MANIFOLD  AND  BOUNDARY  LAYER  STABILIZERS 


Moure  7.7  P I TnT -STATIC  PROBES  AND  TRAVERSE  MECHANISMS 
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Figure  3.8  THE  BUBBLE  GENERATOR  SYSTEM 


Figure  3.10  DISCHARGE  ELBOW  AND  DAMPER 


Release  Tube 


BUBBLE  generator  system  schematic 


3 .4  Calibration  and  Tune-up  of  the  Wind  Tunnel 

Before  making  the  final  test  models,  a  trial  plywood  model  was  fabricated 
and  several  exploratory  runs  were  made  in  the  wind  tunnel  to  ensure  that 

1)  the  required  approach  wind  boundary  layer  velocity  profiles  are 
obtainable  on  the  floor  of  the  test  section, 

2)  boundary  layers  on  the  walls  and  ceiling  of  the  tunnel  are  not 
too  thick, 

3)  with  the  selected  linear  scale  of  1:50,  blockage  effects  due 
to  the  model  are  insignificant,  and 

4)  clearances  between  model  and  tunnel  walls  and  ceiling  are  large 
enough  to  render  errors  from  boundary  layer  interference  effects 
insignificant. 

First,  with  no  model  in  the  test  section,  vertical  traverses  were  made 
with  the  pitot-static  probe  P^  (see  Figures  3.2  and  3.7)  at  span-wise  loca¬ 
tions  from  wall  to  wall  to  determine  velocity  profiles.  The  wall  boundary 
layer  thicknesses  of  the  tunnel  were  observed  to  be  not  more  than  about  3 
inches.  Probe  P2  was  then  moved  along  the  tunnel  axis  to  its  most  forward 
position  and  vertical  traverses  made  with  it.  The  velocity  profiles  were 
identical  to  those  obtained  from  readings  of  probe  P^  in  planes  near  to  the 
tunnel  axis.  Typical  profiles  obtained  from  readings  of  probe  P2  at  a 
counter-jet  manifold  pressure  of  4  inches  of  water  are  shown  in  Figures  3.12, 
3.13,  and  3.14  (for  all  the  test  runs,  counter-jet  manifold  pressure  was 
maintained  at  4  inches  of  water).  It  can  be  seen  that  the  distribution  fits 
very  closely  a  power  law  with  an  exponent  of  1/4.5.  (Recommended  values  of  the 
exponent  for  a  bush  country  vary  between  1/6  and  1/3.5,  Ref.  3.7).  The 
boundary  layer  thickness  is  about  18  inches,  giving  a  model  height  to  boundary 
layer  thickness  ratio  of  1:4.5.  This  translates  to  a  wind  boundary  layer 
thickness  of  72  feet  approaching  a  single  story  building  16  feet  tall  (this 
falls  well  within  the  range  of  60  to  80  feet  recommended  for  open  terrain  or 

bush  country).  Reynolds  numbers  based  on  tunnel  floor  boundary  layer  thick- 

5  5 

ness  range  from  0.7  x  10  to  1.7  x  10  ,  granting  good  turbulent  mixing  in 
these  layers  (the  value  of  critical  Reynolds  number  based  on  boundary  layer 
thickness  for  flow  over  a  flat  plate  is  only  about  3  x  10  ) . 
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Figure  3.13  BOUNDARY  LAYER  VELOCITY  PROFILE  AT  19.3  FPS 


Figure  3.14  BOUNDARY  LAYER  VELOCITY  PROFILE  AT  12.7  FPS 


1 

Next,  the  damper  opening  on  the  discharge  elbow  of  the  blower  was  cali¬ 
brated  against  free  stream  velocity.  As  a  result,  a  single  reading  of  either 
probe  P-|  or  P2  (at  a  height  more  than  18  inches  above  the  tunnel  floor)  was 
sufficient  to  check  the  accuracy  of  the  damper  setting. 

The  trial  model  was  then  placed  inside  the  test  section,  about  20  inches 
downstream  of  probe  P2>  The  overall  dimensions  of  this  model  equalled  those 
of  the  proposed  test  models  (Section  3.5),  which  were  designed  to  a  linear  scale 
of  1:50.  The  blockage  ratio  based  on  projected  frontal  area  was  less  than  4%. 

The  model  heioht  was  4  inches  giving  a  clearance  of  26  inches  between  its  roof 
and  the  tunnel  ceiling.  On  the  sides,  the  clearance  between  the  foot  of  the 
berms  and  the  tunnel  walls  was  20  inches  for  zero-degree  wind  angle. 

With  the  trial  model  located  as  described,  a  vertical  traverse  was  made 
using  probe  P2  which  was  20  inches  upstream  of  the  model.  The  velocity 
profiles  closely  approximated  those  obtained  without  the  model.  Probe  P2  was 
moved  downstream  to  locations  directly  above  the  model  and  both  lateral  and 
vertical  traverses  were  made.  Maximum  velocities  directly  above  the  model 
were  within  5%  of  free  stream  values  showing  that  the  clearance  between  model 
and  tunnel  ceiling  was  large  enough  to  yield  accurate  test  results.  Lateral 
traverses  revealed  that  approach  wind  boundary  layer  velocity  profiles 
obtained  without  the  model  were  approximated  at  locations  between  the  model 
and  tunnel  walls.  This  proved  that  the  side  clearances  were  also  large  enough 
to  yield  accurate  test  results  with  models  of  the  size  selected. 

Several  trial  runs  were  also  made  to  determine  optimum  settings  of  pressure 
regulator  and  metering  valves  of  the  bubble  generator  unit,  locations  and 
maneuverability  of  bubble  release  tubes,  lighting  and  background  settings,  and 
camera  location  and  film  speed.  Extended  trial  runs  were  necessary  before  a 
satisfactory  system  could  be  put  together  that  is  easy  to  operate  and  yields 
acceptable  data  (a  typical  set-up  is  shown  in  Figure  3.9). 

3.5  Description  of  Scaled  Models 

At  the  time  of  our  investigations,  it  was  not  known  what  the  average  upgraded 
fallout  shelter  would  look  like  as  they  had  not  been  located  and  surveyed.  The 
shelters  surveyed  as  part  of  the  National  Fallout  Shelter  Survey  several  years 
ago  did  not  necessarily  represent  the  shelters  now  being  considered.  [CPA 
indicated  that  the  upgraded  shelter  would  probably  be  single-story,  residential 
and  small  commercial  type  structures  that  could  support  full  earth  berms. 
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Since  the  full  scope  of  these  type  structures  could  not  be  completely 
investigated  as  part  of  this  study,  it  was  mutually  decided  to  keep  the 
shelter  configurations  simple  to  ensure  good  accuracy  and  reliability  of  the 
data  taken  during  the  scaled  model  testing.  The  general  configuration  of  a 
bermed  shelter  is  shown  in  Figure  3.15.  The  rectangular  structure  is  typical 
of  a  small  commercial  building,  one  story  high  with  slab-on-grade  construction. 
Openings  can  be  either  windows  or  doors  with  a  variation  in  opening  patterns 
possible.  The  slope  of  the  bermed  sides  depends  on  the  type  of  soil  used, 
varying  from  30°  to  45°.  For  the  scale  models,  an  angle  of  36°  was  used. 

Floor  to  ceiling  wall  partitions  with  doorways  offer  the  main  internal 
impedance  to  airflow  through  the  building. 

As  per  mutual  agreement  it  was  decided  that  the  shelter  configurations 
to  be  studied  should  represent  extreme  conditions,  i.e., 

1)  above  and  below  grade  shelters 

2)  high  and  low  internal  impedance  to  airflow 

3)  simple  and  complex  opening  patterns 

4)  high  and  low  occupant  densities 

Figure  3.16  presents  the  structures  that  contain  this  range  of  character¬ 
istics  and  are  summarized  in  Table  3.1.  Applying  the  principles  of  fluid 
mechanics,  it  is  possible  to  imagine  the  airflow  patterns  that  might  result 
within  each  shelter  configuration  for  various  opening  patterns.  These  are 
depicted  in  Figure  3.17.  For  this  phase  of  the  study,  investigations  were 
limited  to  simple  configurations  with  no  internal  partitions  and  with  scaled 
dimensions  adjusted  to  ensure  no  tunnel  blockage  effects  as  model  was  rotated 
to  simulate  various  wind  angles. 

Figure  3.18  shows  shelter  geometries  modeled  in  the  present  study.  Shelter 
type  I  has  a  door  in  the  front  and  rear,  type  II  additionally  has  a  window 
on  one  of  the  sides  and  type  III  has  windows  on  both  sides.  Models  are  assembled 
from  triangular  pine-wood  prisms  representing  the  earth  berms,  a  plexiglas  sheet 
for  the  floor  and  a  plywood  sheet  for  the  ceiling.  Changing  from  one  model  type 
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Table  3.1 


Summary  of  Shelter  Model  Characteristics 


Shelter 

Configuration 

■Hi 

Opening  Pattern 

Grade 

Location 

Low 

■in 

A1 

✓ 

✓ 

✓ 

A2 

✓ 

/ 

/ 

B1 

/ 

✓ 

/ 

B2 

/ 

/ 

✓ 

Cl 

/ 

/ 

✓ 

C2 

/ 

/ 

/ 

D1 

/ 

/ 

✓ 

D2 

/ 

/ 

✓ 
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to  another  can  be  easily  accomplished  with  this  setup.  The  prisms  and  the  ply¬ 
wood  sheet  are  painted  black  on  all  sides  to  eliminate  undesired  light  reflec¬ 
tions. 

Models  were  made  to  a  linear  scale  of  1:50.  A  full-size  shelter  53  feet 
long  x  40  feet  wide  x  16.5  feet  tall  scaled  to  a  model  12-23/32  inches  x  9-5/8 
inches  x  4  inches.  The  smallest  of  the  model  dimensions  -  namely  the  height 
of  4  inches  -  gives  a  Reynolds  number  of  about  104  at  the  lowest  free-stream 
wind  velocity  tested  (7.5  feet  per  sec.).  This  value  is  well  above  the  upper 
bound  of  5,000  quoted  in  reference  3.1  for  shedding  of  vortices  of  the  Karman 
vortex  street  type. 

3.6  Test  Procedure 

Model  I  was  placed  in  the  test  section  so  that  it  faced  the  approaching 
flow  at  a  prescribed  angle.  Probes  and  P^  were  set  at  20"  above  the  tunnel 

floor.  Two  250-watt  lights  (effective  diameter  about  2")  were  used  to  illumi¬ 
nate  the  bubble  paths.  One  of  the  lights  was  placed  18"  behind  the  rear  door, 
the  other  30"  upstream  of  the  front  door  opening.  The  light  beams  were  directed 
through  the  doors  to  illuminate  the  door-to-door  pathway.  The  damper  opening 
was  set  for  a  particular  free-stream  wind  velocity  and  the  blower  was  turned 
on.  After  conditions  became  steady,  probe  P^  was  read  to  obtain  the  free- 
stream  wind  velocity.  (Probes  P^  and  read  the  same  ensuring  that  blockage 
effects  were  insignificant.)  Thin  white  lines  were  marked  at  1/4"  intervals 
for  a  total  length  of  2  inches  across  each  opening  to  serve  as  distance 
markers. 

A  fine  quality  mirror  (3'  x  2')  mounted  on  an  adjustable  frame  was 
positioned  beneath  the  tunnel  test  section  at  45°  to  the  floor.  The  movie 
camera  was  stationed  at  a  distance  of  about  4  feet  from  the  center  of  the 
mirror,  and  focused  on  the  mirror  image  of  the  model's  interior.  Helium  and 
compressed  air  supplies  to  the  bubble  generator  were  opened  and  the  control 
knobs  of  the  console  adjusted  to  release  neutrally  buoyant  bubbles  ('  1/8- 
inch  diameter)  through  the  release  tubes  at  the  desired  rate.  The  lights 
illuminating  the  doorways  were  turned  on.  Positions  of  the  bubble  release 
tubes  were  adjusted  for  optimal  frequency  of  bubble  flow  through  the  model. 

The  camera  was  turned  on  for  about  30-40  seconds  at  a  nominal  speed  of  64 
frames  per  second.  Thereafter,  the  lights  and  the  bubble  generator  were 
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turned  off  and  the  damper  opening  adjusted  for  a  different  free-stream  wind 
velocity.  The  test  was  repeated  for  the  new  wind  velocity  after  conditions 
became  steady.  For  each  orientation  of  the  building,  bubble  flow  was  filmed 
for  5  different  free-stream  wind  velocities.  Once  this  was  completed,  the 
test  section  access  door  was  opened  and  the  model  rotated  to  a  new  angle  for 
a  different  orientation.  With  each  model,  filming  was  done  for  5  different 
relative  wind  angles  (0°  to  90°). 

When  all  the  measurements  with  model  I  were  completed,  model  II  was 
placed  in  the  test  section.  An  additional  light  source  was  installed  and 
the  beam  from  it  directed  through  the  side  window  opening.  Filming  of  this 
series  was  done  following  the  same  procedure  as  in  the  previous  series. 

After  completing  tests  on  model  II,  model  III  was  placed  in  the  test  section 
and  similar  tests  were  performed.  For  accurate  determination  of  film  speed, 
the  movement  of  the  seconds  needle  of  a  precision  stop  watch  was  also 
periodically  photographed  for  five  seconds.  The  correct  film  speed  was  found 
to  be  73  frames  per  second. 

When  the  planes  of  the  side  openings  were  at  acute  angles  to  the  wind 
direction,  it  was  difficult  to  discern  whether  bubbles  were  entering  or 
leaving  through  the  side  openings.  This  was  due  to  the  complex  vortex  motion 
just  outside  of  these  openings.  In  such  cases,  additional  runs  were  made  with 
the  bubble  release  tubes  turned  sideways.  This  increased  the  external  bubble 
flow  rate  past  the  side  openings  considerably  and  reduced  the  number  of 
bubbles  entering  through  the  front  opening  to  practically  zero.  From  these 
runs  one  could  determine  the  proper  bubble  flow  direction  through  the  side 
openings  without  ambiguity. 

3 ■ 7  Description  of  Re suits 

The  processed  films  were  projected  on  a  screen  frame  by  frame  and  the 
bubble  movement  distance  through  an  opening  per  frame  was  noted.  Knowing  the 
film  speed  (73  frames  per  second)  and  the  distance  moved  per  frame,  the  bubble 
velocity  through  the  opening  could  be  calculated.  The  mean  air  flow  velocity 
through  an  opening  was  obtained  by  averaging  the  velocities  of  15-20  bubbles. 
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Let, 


x  =  bubble  movement  distance  in  inches  through  an  opening  in 
one  frame, 

a  =  cross-sectional  area  of  the  model  opening  in  square  inches, 

n  =  number  of  bubbles  whose  velocity  is  averaged. 

Then,  the  average  air  flow  velocity  (v)  in  feet  per  second 
through  the  model  ODening  is, 

n  xi  1  1 

v  =  X  (1/73)  x  n 

=  6.083  x  x 

n 

•_  xi 

where  v  =  - -  • 

A  n 

The  volume  flow  rate  (a)  in  cubic  feet  per  minute  through  the  model  opening  is, 
q  =  6.083  x  x  x  *  60 
=  2.535  x  x  x  a 

With  a  linear  scale  of  1:50,  flow  rate  through  the  full  size  shelter  opening 
can  be  calculated  as, 

Q  =  q  x  50  x  50 
=  6336.8  x  x  x  a 

When  there  are  more  than  one  inlet  or  exit  openings,  the  throughput  flow  rate 
is  taken  as  the  sum  of  the  flow  rates  through  all  the  inlet  or  exit  openings. 
The  variation  of  ventilation  throughput  with  the  free  stream  wind  velocity  (I4J 
is  shown  in  araphical  form  in  Figures  3.19,  3.20,  and  3.21.  The  same  data  is 
shown  in  tabular  form  in  Tables  3.2,  3.3  and  3.4. 
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Table  3.2 


VENTILATION  THROUGHPUT  CFM  AT 
VARIOUS  WIND  VELOCITIES  AND  ANGLES,  MODEL  I 


1-0° 


WIND  VELOCITY,  FPS 

TOTAL  CFM  -  SCALED 

-  FULL  SIZE 

7.7 

3.5 

3750 

11.3 

4.3 

10750 

14.3 

5.5 

13750 

16.5 

5.8 

14500 

19.1 

6.2 

15500 

1-15° 

WIND  VELOCITY,  FPS 

7.8 

11.6 

14.5 

16.4 

18.8 

TOTAL  CFM  -  SCALED 

3.5 

5.0 

6.2 

6.9 

7.2 

-  FULL  SIZE 

8750 

12500 

15500 

17250 

18000 

1-30° 

WIND  VELOCITY,  FPS 

7.1 

11.2 

14.2 

16.4 

19.1 

TOTAL  CFM  -  SCALED 

3.3 

4.8 

5.6 

6.5 

7.5 

-  FULL  SIZE 

8250 

12000 

14000 

16250 

18750 

1-45° 

WIND  VELOCITY,  FPS 

6.8 

11.0 

14.0 

16.0 

18.8 

TOTAL  DFM  -  SCALED 

3.7 

4.8 

5.9 

6.9 

7.2 

-  FULL  SIZE 

9250 

12000 

14750 

17250 

18000 

1-90° 

VIRTUALLY  NO  FLOW  AT  ALL  WIND  VELOCITIES. 
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Table  3.3 


VENTILATION  THROUGHPUT  CFM  AT 
VARIOUS  WIND  VELOCITIES  AND  ANGLES,  MODEL  II 


11-0° 

WIND  VELOCITY,  FPS 

TOTAL  CFM  -  SCALED 

-  FULL  SIZE 

7.3 

3.9 

9750 

11.5 

5.2 

13000 

13.7 

6.9 

17250 

15.9 

7.1 

17750 

19.0 

8.2 

20500 

11-15° 

WIND  VELOCITY,  FPS 

7.8 

11.6 

13.9 

15.8 

19.1 

TOTAL  CFM  -  SCALED 

4.6 

5.7 

6.0 

7.4 

8.4 

-  FULL  SIZE 

11500 

14250 

15000 

18500 

21000 

11-30° 

WIND  VELOCITY,  FPS 

7.4 

11.5 

14.4 

16.2 

18.7 

TOTAL  CFM  -  SCALED 

3.9 

5.4 

6.3 

7.4 

7.8 

-  FULL  SIZE 

9750 

13500 

15750 

18500 

19500 

►—4 

1 

cn 

o 

WIND  VELOCITY,  FPS 

7.1 

11.6 

13.9 

16.2 

18.7 

TOTAL  CFM  -  SCALED 

4.4 

5.5 

6.8 

7.3 

8.7  i 

-  FULL  SIZE 

11000 

13750 

17000 

18250 

21750 

11-90° 

WIND  VELOCITY,  FPS 

7.3 

12.0 

14.2 

16.3 

19.2 

TOTAL  CFM  -  SCALED 

6.4 

8.7 

10.5 

10.8 

12.0 

-  FULL  SIZE 

16000 

21750 

26250 

27000 

30000 
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Table  3.4 


VENTILATION  THROUGHPUT  CFM  AT 
VARIOUS  WIND  VELOCITIES  AND  ANGLES,  MODEL  III 


111-0° 

WIND  VELOCITY,  FPS 

TOTAL  CFM  -  SCALED 

-  FULL  SIZE 

7.3 

4.3 

10750 

11.6 

6.5 

16250 

14.5 

7.2 

18000 

16.5 

7.7 

19250 

18.6 

8.7 

21750 

111-30° 

WIND  VELOCITY,  FPS 

8.1 

12.7 

15.7 

18.7 

TOTAL  CFM  -  SCALED 

7.5 

11.1 

14.7 

16.6 

-  FULL  SIZE 

18750 

27750 

36750 

41cCo 

111-45° 

WIND  VELOCITY,  FPS 

7.3 

11.7 

14.6 

16.2 

18.7 

TOTAL  CFM  -  SCALED 

8.8 

13.1 

16.2 

17.5 

19.0 

-  FULL  SIZE 

22000 

32750 

40500 

43750 

47500 

111-90° 

WIND  VELOCITY,  FPS 

6.6 

11. 1 

14.1 

16.0 

18.8 

TOTAL  CFM  -  SCALED 

8.7 

9.4 

10.4 

12.2 

13.6 

-  FULL  SIZE 

21750 

23500 

26000 

30500 

34000 
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3 . 8  Discussion  of  Results 


Several  important  observations  and  conclusions  can  be  made  from  this  study. 

These  are  discussed  below. 

1)  At  any  given  relative  wind  angle  e,  the  ventilation  flow  rate  (q) 
through  all  three  models  varies  linearly  with  wind  speed  (see 
Figures  3.1b,  3.20  and  3.21). 

2)  For  models  I  and  II,  the  ventilation  flow  rate  does  not  vary 
significantly  with  relative  wind  angle  for  values  of  e  not  too 
large  (as  can  be  seen  from  the  results  for  0°  <  e  <45°).  Flows 
at  non-zero  values  of  e,  tend  to  be  slightly  larger  than  those 
at  e  =  0°,  because  of  the  higher  average  pressure  in  the  stagna¬ 
tion  zone  immediately  upstream  of  the  front  door  created  by  the 
sides  of  the  berm.  But  this  tendency  is  counteracted  by  the 
smaller  door  area  normal  to  the  wind  at  non-zero  angles.  (For 
model  I,  the  total  normal  area  decreases  with  wind  angle.  For 
model  II,  the  opening  area  normal  to  the  wind  steadily  increases 
as  0  goes  from  0°  to  90°).  For  small  values  of  e,  the  former 
factor  dominates  and  for  larger  values  of  0  the  latter  factor  is 
dominant.  A  comparison  of  model  I  results  for  e  _<45°  (Figure 
3.19)  against  those  for  90°  (and  80°)  would  reveal  this.  Although 
filming  was  not  done,  bubble  flow  through  model  I  was  observed  at 
e  =  80°.  A  very  small,  but  significant,  through  flow  was  noticed. 

For  model  III,  throughput  flow  rate  at  non-zero  values  of  e  is 
considerably  greater  than  that  for  e  =  0°.  This  is  what  one  would 
predict  noting  that  both  the  average  stagnation  pressure  and  the 
total  projected  normal  opening  area  are  larger  for  winds  at  non¬ 
zero  angles.  At  8  =  45°,  for  which  the  projected  normal  area  is 
maximum,  throughput  flow  rate  for  model  III  has  its  largest  value. 

3)  Increasing  the  opening  area  on  the  windward  side  alone  or  on  the 
leeward  side  alone  does  not  augment  the  throughput  flow  rate  very 
much.  To  realize  maximum  throughput  flow  rate,  the  total  opening 
area  must  be  more  or  less  equally  divided  between  the  windward  and 
leeward  sides.  A  comparison  of  the  results  for  models  I,  II  and 
III  clearly  bears  this  out. 
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Model  II  has  about  74%  more  total  opening  area  than  model  I  and 
model  III  has  43%  more  opening  area  than  model  II.  For  e  =  15°, 
30°  and  45°,  throughput  flow  rate  for  model  II  is  only  about  20% 
larger  than  for  model  I,  whereas  for  model  III  it  is  about  80% 
more  than  that  for  model  II.  Although  model  II  has  considerably 
larger  area  than  model  I  on  the  windward  side,  it  has  the  same 
area  as  model  I  on  the  leeward  side.  This  severe  flow  restriction 
on  the  leeward  side  prevents  the  model  II  flow  from  going  sub¬ 
stantially  higher  than  that  for  model  I.  Model  III,  on  the  other 
hand,  has  an  even  distribution  of  opening  areas  between  the  wind¬ 
ward  and  leeward  sides.  It  is  also  of  interest  to  note  that, 
model  II  has  a  much  larger  projected  area  of  the  inlet  openings 
normal  to  the  wind  at  e  =  45°;  yet,  the  throughput  flow  rate  is 
considerably  smaller  than  for  e  =  90°.  This  is  due  to  the  better 
distribution  of  opening  areas  between  windward  and  leeward  sides 
in  the  latter  case.  As  would  be  expected,  the  highest  flow  rate 
for  model  III  occurs  at  a  relative  wind  angle  of  45°. 

4)  For  model  I  at  90°  relative  wind  angle,  there  is  practically  zero 
flow  through  the  openings  (for  this  case  both  door  openings  are 
in  planes  parallel  to  the  wind  velocity).  Model  II,  at  the  same 
wind  angle,  has  some  flow  through  the  side  openings  induced  by 
the  suction  at  the  rear  window  opening. 

5)  Most  of  the  flow  through  model  III  at  90°  relative  wind  angle, 
enters  through  the  front  window  and  exits  through  the  rear  window. 
The  flow  pattern  at  the  side  door  openings  is  very  complex. 

Through  part  of  the  side  openings,  flow  is  going  out  and  through 
some  other  parts  flow  is  coming  in.  The  net  flow  is  outwards 

and  is  calculated  as  the  difference  in  volume  flow  rates  through 
the  front  and  rear  windows. 

It  may  be  noted  that  at  e  =  90°,  model  III  is  essentially  the 
same  as  it  is  at  0°  with  the  only  difference  that  the  opening 
area  normal  to  the  wind  is  greater  by  a  factor  of  1.48.  The 
projected  flow  rates  obtained  by  multiplying  those  at  0°  by 
1.48  are  shown  by  the  chain  line  in  Figure  3.21.  The  measured 
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values  are  about  10%  higher  than  the  projected  values.  This  is 
not  surprising  since  the  length  of  the  main  flow  path  (window-to- 
window  distance)  in  this  case  is  about  30%  less  than  that  for  the 
0°  position  (door-to-door  distance). 

6)  Figures  3.22,  3.23,  3.24,  and  3.25  show  a  comparison  of  inflow 
rates  through  the  front  door  opening  alone  for  models  I,  II  and 
III  for  different  wind  speeds  and  angles.  Several  interesting 
observations  can  be  made  from  these  figures.  At  e  =  0°,  flow 
rates  through  model  III  are  larger  than  those  through  model  II 
which  in  turn  are  slightly  larger  than  flow  rates  through  model 
I.  This  is  to  be  expected  because  of  the  small,  but  non-zero, 
outflow  through  the  side  openings  of  models  II  and  III.  At 
e  =  15°,  there  is  some  inflow  through  the  window  of  model  II. 

The  net  flow  through  the  window  is  still  outwards.  Therefore, 
inflow  through  the  front  door  is  still  higher  than  that  for 
model  I.  As  e  is  increased  to  30°,  the  flow  through  model  II 
window  becomes  unidirectional  and  inwards.  Model  II  for  this 
angle  has  a  much  larger  inlet  opening  area  although  it  has  the 
same  exit  opening  area  as  model  I.  The  increased  flow  causes 
a  pressure  build-up  upstream  of  the  rear  door  inside  the  model. 
Consequently,  the  inlet  flow  through  the  front  door  opening  sees 
a  larger  resistance  and  drops  to  a  value  less  than  that  for  model 
I,  although  the  total  throughput  flow  rate  for  model  II  is  larger 
(see  Figures  3.19  and  3,20).  In  the  case  of  model  III,  this  sort 
of  a  pressure  build-up  does  not  occur  since  the  leeward  side  and 
windward  side  opening  areas  always  match.  However,  as  the  cross 
stream  from  the  windward  side  window  becomes  stronger,  the 
pressure  in  the  mixing  zone  of  the  two  streams  begins  to  rise. 
Consequently,  the  inlet  flow  through  the  front  door  sees  an 
increasing  resistance  and  drops  to  values  slightly  below  those 
for  model  I  (see  results  for  e  =  45°). 
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FULL  SIZE  DOOR  AIRFLO'I  RATE  -  CFM 


SCALED  MODEL  AIRFLOW 


r inure  T.DT  I  Nil  ON  RATES  THROUGH  THE  FRONT  DOOR  AT  15° 
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FULL  SIZE  DOOR  AIRFLOW  RATE  -  CFM 


SCALED  MODEL  flIRFLO"  THROUGH  FRO'IT  DOOR  -  CFM 
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CALEU  MODEL  A  I RFLO'.'  THPQMGH  FRONT  DOOR 


FULL  SIZE  DOOR  AIRFLOW  RATE  -  CFM 


Measurement  of  ventilation  throughput  in  model  buildings  poses  several 
difficulties.  Due  to  practical  considerations  the  maximum  dimensions  of  the 
model  and  its  openings  are  limited.  Further,  the  flow  patterns  inside  the 
model  and  in  the  neighborhood  of  the  openings  are  so  complex  that  measurement 
of  flow  through  these  openings  using  probes  is  undesirable.  The  method  of  flow 
tracing  with  neutrally  buoyant  bubbles  and  determining  their  velocity  with 
motion  photography  developed  by  GARD  and  used  in  this  study,  is  effective  and 
relatively  simple.  Measurements  are  repeatable  within  +  5%.  However,  there 
is  one  major  drawback  for  this  method  (GARD  is  already  in  the  process  of  finding 
correction  factors  to  reduce  errors  from  this  source)  which  relates  to  the 
assumption  that  the  average  flow  velocity  obtained  is  a  true  average  for  the 
whole  area  of  the  opening.  In  reducing  data  from  the  films,  for  some  of  the 
runs,  the  average  velocities  of  15-20  bubbles  were  calculated  from  two  different 
segments  of  the  film.  The  maximum  difference  between  the  average  velocities  of 
the  two  sets  was  less  than  10%.  This  shows  that  the  average  velocity  obtained 
is  the  true  average  for  the  "through-flow  area"  of  the  opening.  However,  this 
need  not  be  the  average  velocity  for  the  whole  area  of  the  opening.  A  part  of 
the  opening  area  may  be  effectively  blocked  by  recirculating  vortices,  so  that 
the  actual  through-flow  area  of  the  opening  could  be  significantly  less  than 
its  full  cross-sectional  area. 

An  approximate  value  of  the  "through-flow  area"  for  the  front  door  opening 
of  model  I  can  be  calculated  in  the  following  way.  Let 

Ac  =  coefficient  of  area 

_  through-flow  area  of  opening 
cross-sectional  area  of  opening 

Considering  the  entrance  flow  geometry,  it  appears  reasonable  to  assume  that 
the  coefficient  of  area  of  the  rear  door  opening  ~  unity.  Based  on  this  and 
using  the  law  of  mass  conservation,  values  of  A£  for  the  front  door  opening 
can  be  calculated  as 

A  _  average  flow  velocity  through  rear  door _ 

c  average  velocity  for  through-flow  area  of  front  door 

average  bubble  flow  velocity  through  rear  door 
average  bubble  flow  velocity  through  front  door 
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Values  of  Ac  as  a  function  of  wind  speed  and  angle  are  shown  in  Figure  3.26. 
For  a  given  wind  speed  and  angle,  the  value  of  Ac  for  models  II  and  III  would 
be  nearly  the  same  as  that  for  model  I.  Any  difference  in  the  values  of  Ac 
among  the  different  models  should  arise  from  a  change  in  the  distribution  of 
the  pressure  differentials  across  the  front  door,  resulting  from  the  air  flow 
through  the  windows.  This  effect  is  likely  to  become  significant  only  at 
larger  values  of  o  when  the  momentum  and  mass  flux  through  the  window  openings 
become  large.  Values  of  flow  rate  for  the  front  door  openings  of  models  II 
and  III  plotted  in  Figures  3.22  through  3.25  were  obtained  by  multiplying  the 
product  of  the  measured  average  velocity  and  the  door  area  by  the  coefficient 
of  area  obtained  for  model  I. 

An  important  fact  brought  out  by  the  present  study  relates  to  the  non- 
linearities  inherent  in  fluid  dynamics  problems  of  this  sort  and  the  need  to 
exercise  caution  in  interpolating/extrapolating  experimental  data.  At  e  =  30° 
and  45°,  flow  through  the  side  window  of  model  II  is  inwards.  Extrapolating, 
one  finds  that  the  flow  through  the  window  at  e  =  15°  should  be  smaller,  but, 
still  inwards.  In  reality,  the  net  flow  through  the  model  II  window  at 
o  =  15°  is  outwards  although  there  is  a  very  small  inflow  through  a  section 
of  this  window. 

Another  illustration  of  the  above  facts  can  be  seen  by  attempting  to  pre¬ 
dict  the  throughput  flow  rate  for  model  III  at  e  =  90°,  from  the  results  of 
models  I  and  II.  For  model  III  at  e  =  90°,  one  may  be  tempted  to  calculate 
the  influx  through  the  windward  window  opening  by  multiplying  the  throughput 
for  model  I  at  0°  by  the  area  ratio  of  window  to  door.  Adding  to  this  the 
influx  through  the  side  door  openings  of  model  II  at  e  =  -90°  (although  not 
filmed,  flow  through  model  II  for  e  =  -90°  was  visually  observed,  a  small 
but  significant  flow  enters  throuoh  the  side  doors  and  exits  throuoh  the  rear 
window),  one  might  predict  that  for  model  III  at  e  =  90°,  the  efflux  through 
the  leeward  window  is  more  than  the  influx  through  the  windward  window  and  the 
net  flow  through  the  side  door  openings  is  inwards.  This  is  incorrect  based 
on  visual  observations  and  the  measured  values  of  throughput  flow  rate  for 
model  III  at  o  =  90°. 
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Figure  h26  COEFFICIENT  OF  AREA,  A  ,  VS.  V. 
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Based  on  the  test  results  obtained  from  this  program,  it  is  clear  that 
only  an  extended  scaled  model  study  yielding  correlations  of  throughput  flow 
with  the  dominant  independent  variables  can  provide  the  information  required 
to  make  reasonably  accurate  predictions  and  meaningful  recommendations  with 
respect  to  ventilation  in  shelter  buildings.  This  study  must  include  all  the 
dominant  independent  variables  and  extend  over  a  wide  enough  range  so  that  the 
results  can  be  extrapolated  or  interpolated  with  confidence.  The  dominant 
variables  for  this  study  can  be  identified  as  wind  speed  and  direction;  height, 
width  and  length  of  the  building;  areas  and  locations  of  openings  and  internal 
resistances  including  those  due  to  occupants.  Since  the  number  of  independent 
variables  is  large  and  their  interactions  are  non-linear,  the  testmatrix.  will 
have  to  be  selected  judiciously  to  render  the  study  sufficiently  thorough  and 
detailed  while  practical  to  perform. 
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Section  4 


NATURAL  VENTILATION  ADEQUACY  ASSESSMENT  TECHNIQUE 


To  accomplish  the  ultimate  goal  of  this  work  unit,  i.e.,  determine  the 
adequacy  of  natural  ventilation  for  upgraded  shelters,  the  following  approach 
was  adopted  as  a  means  of  performing  the  adequacy  assessment  calculations: 

1)  Determine  experimentally  the  ventilation  airflow  characteristics 
of  upgraded  shelters  as  a  function  of  such  variables  as  ambient 
wind  speed  and  direction,  shelter  configuration  and  opening 
pattern,  internal  flow  impedance,  etc. 

2)  Develop  a  mathematical  model,  based  upon  the  results  of  (1)  which 
can  be  combined  with  algorithms  for  performing  shelter  heat  balances 
to  determine  resultant  daily  average  shelter  effective  temperatures 
on  the  basis  of  historical  weather. 

3)  Computerize  the  adequacy  assessment  methodology  to  enable  parametric 
studies  for  various  types  of  shelters  and  climatic  zones  to  be  done 
to  determine  if  the  82°F  effective  temperature  and  90%  adequacy 
criteria  can  be  met  for  a  range  of  occupant  densities. 

Step  1,  as  described  in  Section  3,  was  performed  for  a  limited  sample  of 
simple  shelter  configurations  and  limited  range  of  wind  speeds  and  directions. 
Steps  2  and  3  are  described  in  subsequent  portions  of  this  section  along  with 
other  concepts  that  form  the  basis  for  the  adequacy  assessment  technique. 

4. 1  Adequacy  and  Effective  Temperature 

"Effective  temperature"  (TEF)  is  an  empirical  thermal  index  which  combines 
the  effects  of  dry-bulb  temperature,  wet-bulb  temperature  and  air  movement 
to  yield  equal  sensations  of  warmth  or  cold  to  a  person.  "Adequacy"  may  be 
defined  as  the  percentage  of  time  over  some  given  time  domain,  that  the  daily 
average  effective  temperature  of  a  shelter  can  be  expected  to  be  less  than  a 
prescribed  limit.  The  relation  between  the  two  can  perhaps  be  better  seen 
through  an  example. 

Let  us  assume  that  the  shelter  configuration,  the  hourly  weather  data  and 
the  density  of  occupants  are  known,  and  that  the  resulting  average  daily 
effective  temperature  (TEF)  inside  the  shelter  due  to  varying  ambient  weather 
can  be  plotted  against  time  as  shown  in  Figure  4.1.  Also,  let  TEF^  TEF^  and 
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TEF3  be  the  prescribed  limits  for  which  it  is  desired  to  determine  the  per¬ 
centage  of  time  that  the  shelter  was  adequately  ventilated.  For  each  of  these 
temperature  limits,  the  adequacy  (%)  is  the  ratio  of  the  number  of  days  the 
daily  average  inside  effective  temperature  TEF  is  less  than  the  prescribed 
limit.  The  diagram  indicates  therefore  that: 

1)  For  TEF.J,  the  adequacy  is  zero,  since  the  daily  average  effective 
temperature  is  higher  than  TEF^  during  the  entire  period  from  to 


2)  For  TEF^,  the  adequacy,  ADEQ,  is  equal  to 


ADEQ  = 


1  - 


(T,  -  T-i)  +  (T, 


V 


Tu 


V 


x  100 


(Eq.  4.1) 


3)  For  TEF3,  the  adequacy  is  equal  to  100%  since  the  daily  average 
effective  temperatures  are  always  lower  than  the  prescribed  limit. 

Adequacy,  then,  may  be  calculated  for  a  series  of  effective  temperature 
limits.  The  result  is  a  relationship  like  that  depicted  in  Figure  4.2. 

There  will  always  be  effective  temperature  limits,  which  if  set  high  enough 
(TEFb)  or  low  enough  (TEFg),  will  result  in  adequacies  of  100%  and  0%  respec¬ 
tively.  For  TEF^  <  TEF  <  TEF^,  the  adequacy  will  range  between  zero  and  100%. 

The  shape  of  the  curve  and  the  temperature  range  over  which  this  occurs  will 
be  a  function  of  climatic  zone,  shelter  characteristics,  occupant  densities,  etc. 

If  the  effects  of  various  occupant  densities  and  effective  temperature 
limits  on  adequacy  are  also  to  be  considered,  a  family  of  curves  similar  to 
those  shown  in  Figure  4.3  will  result.  The  curve  labeled  n.j  would  represent  the 
unoccupied  condition  while  the  curve  labeled  n^  would  represent  the  maximum 
occupant  density,  which  for  the  purpose  of  this  study  is  taken  to  be  one  person 
per  5  square  feet  or  0.2  person  per  sq.  ft.  The  existance  of  a  relationship 
like  that  depicted  in  Figure  4.3  is  the  basis  of  the  adequacy  assessment  technique,  j 

The  effective  temperature  which  results  inside  a  shelter  depends  on  many  \ 

factors,  includinq: 

•  ventilation  air  flow  rate  (L)  through  the  shelter, 

a  dry-bulb  temperature  (DBT)  and  wet-bulb  temperature  (WBT)  of 
ambient  air, 

t  metabolic  heat  generation  rate  of  occupants  (Qocc)» 

•  heat  transferred  through  the  envelope  of  the  structure, 
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ADEQUACY  (%) 


TER,  TER. 

a  b 

DAILY  AVERAGE  EFFECTIVE  TEMPERATURE  LIMIT 


) 


Figure  4.2  TYPICAL  RELATIONSHIP  BETWEEN'  ADEQUACY  AND 
SHELTER  DAILY  AVERAGE  EFFECTIVE  TEMPERATURE 

LIMIT 
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ADEQUACY  (%) 


EFFECTIVE  TEMPERATURE  LIMIT 


Figure  4.3  TYPICAL  RELATIONSHIP  BETWEEN  OCCUPANT  DENSITY 
AND  EFFECTIVE  TEMPERATURE  LIMIT  ON  ADEQUACY 
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•  additional  heat  loads  like  those  from  direct  sun  light,  electric 
1 i oh t  hulbs  ,  etc. ,  and 

•  air  flow  pattern  inside  the  shelter. 

Little  is  known  presently  about  the  construction  characteristics  of  upgraded 
shelters,  and  therefore  heat  transferred  through  shelter  boundaries  is  assumed 
steady  state  and  additional  internal  heat  loads  other  than  people,  are  being 
ignored.  Further,  it  is  assumed  that  air  within  the  shelter  is  uniformly 
distributed  such  that  no  air  stratification  or  temperature  variation  occurs. 
The  shelter  effective  temperature  is  then  a  function  of  the  following: 

TEF  =  f(()  L,  OPT,  WPT }  (Eq.  4.2) 

Ovt  « 

4 . 2  Occupant  Metabol  i_c  Heat 

The  metabolic  heat  produced  by  an  individual  depends  on  the  sex,  age  and 
activity  level.  Average  metabolic  rates  for  sedentary  activity  levels  as 
functions  of  age  and  sex  for  the  most  recent  census  of  U.S.  population  are 
given  in  the  1 97S  ASHRAE  Applications  Handbook  (Ref  4.1,  Fig.  1.  page  12.2). 
Table  4.1  below  was  developed  using  these  data. 

Table  4.1 

Sedentarv  Metabolic  Rates  as  a  Function  of  Age  and  Sex 


j  PERCENT 

Or  POPULATION*  I 

METALLIC 

HEAT  BTUH 

AGE 

P  MALE 

i  female  i 

MALE 

FEMALE 

75  and 

over  1.6 

310 

200 

70  to 

74 

1  . 1 

1.5 

320 

210 

55  to 

60 

1.5 

1  .8 

350 

230 

60  to 

64 

■> 

3.3  ! 

360 

235 

55  to 

60 

3 . 3 

3.7 

370 

240 

50  to 

54 

3.6 

3 . 8 

375 

245 

45  to 

40 

3  0 

3.1 

300 

250 

40  to 

44 

2 . 0 

3 

30  6 

250 

T6  to 

.10 

>  7 

1 

3.0  i 

400 

250 

30  to 

34 

•  ^  0 

1  *-  • 1 

! 

400 

250 

26  to 

30 

^ 

3.3 

400 

250 

20  to 

34 

3.0 

4.3  : 

405 

260 

16  to 

10 

4.7 

4.6  j 

400 

260 

10  to 

14 

5.3 

5 

375 

250 

5  to 

0 

5 

4.8 

300 

200 

Under 

5 

4.3 

4.1 

1 50  ; 

i 

150 

Total 

p  43.0 

51.3 

l  ...  . 

I 

1 

- 

*  "Social 

Indicators," 

1073,  U.  S.  Uept. 

of  Commerce 
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Houghten  (Ref.  4.2)  derived  polynomial  expressions  for  metabolic  heat 
generation  that  are  convenient  for  computerized  calculations  and  are  as  follows: 


^occ,  total 
^occ,  sensible 


a(TEF;  +  b(TEF)  +  c 
A(TEF)2  +  B(TEF)  +  C 


(Eq.  4.3) 


where 


1)  a  -  0;  b  =  -1.482;  c  =  514.  for  50°  -  TEF  <  87°F 

2)  a  =  -1.508;  b  =  259.7;  c  =  -10795.2.  for  87°  v  TEF  <  102°F 

3)  a  =  0;  b  =  0;  c  *  0,  for  TEF  -  102°F 

4)  A  =  -0.06875;  B  =  1.625;  C  -  523.  for  DBT  "  50°F 


occ,  total 


0  +  n 

■occ,  sensible  occ.  latent 


(Eq.  4.4) 


The  applicability  of  these  equations  have  been  proved  for  only  healthy  male 
subjects  in  the  aqe  group  18-24  years.  For  this  study  an  "average  person" 
whose  metabolic  heat  rate  is  the  weighted  average  of  that  given  in  Table  4.1 
has  been  created.  Osina  the  following  relation: 


'average  person 


('■  of  males)  x  (BTOH  male)  + 

=  of  females)  x  (BTOH  female)  _  p 


(Eq.  4.5) 


This  is  70''-  of  the  value  given  in  the  table,  for  males  between  20-24  years. 

Its  inn  these  equations  and  data  from  ASH  RAF  handbook  of  Fundamentals  (Ref.  4.3, 
page  8-6,  Tables  2  and  3),  the  following  equations  were  generated  for  cal¬ 
culation  occupant  metabolic  heat: 


occ,  total 


.413  a(TEF)  +  b(TEF)  +  c 


1 

N  + 


1.091  >’  V.u 


occ,  sensible 


=  .413 


A(TFF)1-  +  B(TEF)  +  C  N  +  1.091  r.  V.( 


Eq.  4.6) 


where  N 


number  of  occupants 

air  velocity  felt  by  occupant  -  i 

effective  temperature  of  shelter  air  and  can  be  calculated  as 


107. 5(1 OPTS)  -  45. 2(WRTS) 
TOOTS)'  -  t'WBTST  +  62.3 


(Eq.  4.7) 


where  OPTS  -  dry-bulb  temperature  of  shelter  air 
WBT5  -  wet-bulb  temperature  of  shelter  air. 


GARD.  INC. 


4-7 


4.3  Shelter  Heat  Balance 

A  fallout  shelter  is  a  volume  inclosed  by  boundaries  into  which  sensible 
and  latent  heat  loads  and  ventilating  air  are  introduced  and  from  which  air 
is  exhausted  and  energy  is  lost,  see  Figure  4.4.  The  following  assumptions 
are  made: 

1)  the  air  within  the  shelter  is  completely  and  instantaneously  mixed, 

2)  the  film  heat  transfer  coefficients  are  constant  for  any  one 
boundary  surface, 

3)  the  radiative  energy  transfer  within  the  shelter  can  be  neglected, 

4)  the  condition  of  the  air  exhausted  from  the  shelter  is  the  con¬ 
dition  of  the  shelter  atmosphere, 

5)  the  thermal  and  physical  properties  of  the  structural  materials 
and  of  air  are  not  temperature-dependent, 

6)  the  incident  solar  radiation  is  absorbed  on  the  outer  surface  of  the 
boundaries  and  appears  as  conducted  energy  or  is  transmitted  into 
the  shelter  and  is  considered  an  instantaneous  load  along  with  other 
internal  loads,  and 

7)  the  thermal  loads  and  the  psychrometri c  states  of  the  inlet  and 
exhaust  air  are  constant  over  short  time  intervals. 

A  heat  balance  of  the  shelter  volume  yields  the  following: 

SENSIBLE  HEAT  GAIN  =  Qocc>  sensible  +  K  x  F  x  (DBT  -  DBTS) 

LATENT  HEAT  GAIN  =  f )  ,  .. 

occ,  latent 

SENSIBLE  HEAT  LOSS  =  60xc  xpxLx  (DBTS  -  DBT)  ^Eq'  4'8^ 

LATENT  HEAT  LOSS  =  0.625  x  10'5  x  r  x  L  x  (WS  -  W) 


Since  the  shelter  must  be  in  thermal  equilbirium  at  the  end  of  each  time 
interval , 

HEAT  GAIN  =  HEAT  LOSS 


Therefore, 

Q 


and  Q 


=  (60  xcxpxL  +  Kx  F)(DBTS  -  DBT) 
occ,  sensible  '  p 

1-5 


=  0.625  x  10'J  x  p  x  L  (WS  -  W) 


(Eq.  4.9) 


occ,  latent 
where  p  =  density  of  leaving  air 
L  =  ventilation  throughput 
K  =  overall  shelter  heat  transfer  coefficient 
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Figure  4.4  TYPICAL  HEAT  BALANCE 
REQUIRED  FOR  SHELTERS 
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F  =  total  external  surface  area  of  shelter 
W$  =  humidity  ratio  of  shelter  air 
W  =  humidity  ratio  of  outside  air 
cp  =  specific  heat  of  air 
Using  equations  4.3  through  4.9,  one  can  obtain 


-(N  x  B  +  Z) 


[H  x  B  -  Z)  -  4(N  x  A)  Z(DBT)  +  N  x  C 
2  N  x  A 


(Eg.  4.10) 


where  Z  =  0.000355  x  (PB)  x  L  +  K  x  F 

N 

N  =  0.413  (N  +  1.091  E  V. ,D) 

i  =  l  1 

and  PB  =  barometric  pressure 

Once  the  ventilation  throughput  rate,  L,  is  established  for  each  time  interval, 
the  resulting  shelter  effective  temperature  can  be  calculated  for  a  range  of 
occupant  densities  using  the  above  equation. 

4.4  Ventilation  Air  Throughput 

Two  different  approaches  can  be  taken  in  setting  up  the  ventilation  air 
throughput  calculation  model,  i.e.,  (1)  a  global  approach  yielding  a  single 
general  equation  and  (2)  a  differential  approach  yielding  a  system  of  equations, 
one  for  each  of  the  openings.  The  former  is  preferred  and  results  in  a  single 
relationship  which  can  be  solved  once  empirical  constants  have  been  determined 
by  scaled  model  testing. 

The  variables  that  affect  the  ventilation  air  throughput  rate  can  be 
identified  as, 

•  building  configuration, 

•  location  of  openings, 

•  wind  speed  and  direction, 

•  area  of  openings, 

•  ratio  of  openings  on  windward  and  leeward  sides,  and 

•  shelter  orientation. 

Stack  effect  is  ignored  since  it  is  likely  to  be  important  only  in  tall  buildings 
or  under  no  wind  conditions. 

The  ventilation  air  flow  rate  (L)  may  be  expressed  as  (see  Appendix  B 
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(Eq.  4.11) 


for  theoretical  derivation): 

L  =  m(sc,e,v)  xSnFex 

where 


m(sc  0  v)  =  shelter  air  resistance  coefficient 

S  =  specific  area  coefficient 
£Fe  =  sum  of  all  external  wall  opening  areas 


g  =  gravitational  acceleration 

h  =  "pressure  head"  difference  between  the  windward  and  leeward 
sides  of  the  shelter  at  a  "standard"  wind  velocity  and 
direction  for  a  particular  shelter  configuration  and  is 
assumed  to  be  0.5  ft.  of  air 


The  specific  area  coefficient,  S,  is  defined  as. 


J  =  2,4,6, ...,2N 


aj  =  zfe 

Aj  =  area  of  jtFl  opening 


(Eq.  4.12) 


where 

N  =  Number  of  flow  resistances  between  the  windward  and  leeward 
sides. 

The  function  m(SCj0tV)  is  a  function  of  the  shelter  configuration 
(scheme)  -  sc,  the  relative  wind  angle  between  the  wind  and  the  shelter 
orientation  -  0,  and  the  wind  velocity  -v.  The  empirical  function,  m,  can 
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only  be  determined  by  means  of  scaled  model  testing  due  to  its  complex  nature. 
To  accomplish  this,  equation,  4.11  was  rearranged  as  follows. 

„  L  _ 

(Eq.  4.13) 


sc,e,v 


S  x  (EFe)  x  -y2g  x  h  ^ 


Scaled  model  testing  as  described  in  Section  3  was  performed  to  determine  the 
ventilation  air  throughput,  L,  as  a  function  of  the  shelter  scheme,  sc, 
relative  wind  angle,  0,  and  wind  velocity,  v.  Coefficients  S  and  EFe  were 
calculated  from  the  shelter  geometry.  The  pressure  head  difference,  h,  was 
assumed  to  be  equivalent  to  0.5  ft.  air.  The  empirical  function,  m,  was  then 
calculated. 

Experimental  data  was  obtained  for  three  shelter  schemes  at  relative 
wind  angles  of  0°,  15°,  30°,  45°  and  90°  and  for  wind  velocities  of  6-19  feet 
per  second  (4-13  mph).  The  resultant  values  of  m  were  thus  calculated  and 
incorporated  into  the  computerized  adequacy  assessment  model. 

4.5  Adequacy  Assessment  Technique 

A  computerized  adequacy  assessment  technique  was  developed  to  allow 
quantitative  determination  of  adequacy  for  a  given  shelter  and  occupant  loading 
in  a  given  climatic  zone  on  the  basis  of  historical  weather.  A  computerized 
methodology  also  gives  the  capability  of  easily  performing  parametric  analysis. 
A  flowchart  of  the  assessment  methodology  is  shown  in  Figure  4.5.  For  a 
specified  shelter  scheme,  shelter  orientation,  occupant  density,  and  climatic 
zone,  weather  data  is  read  periodically  from  a  historical  file  and  the  shelter 
ventilation  air  throughput  rate  calculated.  A  total  shelter  heat  balance  is 
then  performed  and  the  resultant  shelter  effective  temperature  determined. 

The  resultant  daily  average  effective  temperature  is  tabulated  over  the  period 
of  analysis  and  the  resultant  adequacy  calculated.  The  process  can  be  repeated 
for  other  densities,  orientations  and  shelter  schemes. 

Historical  weather  data  from  the  National  Oceanic  and  Atmospheric  Ad¬ 
ministration  is  used  as  a  source  of  hourly  data  for: 

1)  dry-bulb  temperature 

2)  wet-bulL  temperature 

3)  wind  velocity 

4)  wind  direction 

5)  barometrc  pressure 

6)  cloud  cover 


WEATHER 


WEATHER  DATA  DECODING 


WEATHER 

TAPE 


As  per  FEMA's  recommendation,  the  period  of  analysis  was  limited  to  the  10% 
hottest  days  of  the  year,  which  it  was  assumed  fell  within  the  period  from 
May  1  throuqh  September  30th.  Samples  of  program  output  reports  are  shown  in 
Figures  4.6  to  4. 10. 

To  demonstrate  the  methodology  an  analysis  was  performed  for  Model  I  des¬ 
cribed  in  Section  3,  with  a  southern  orientation  and  two  different  climatic 
zones  represented  by  Atlanta,  Georgia  and  Portland,  Oregon.  The  results  were 
as  follows: 


Table  4.2 


OCCUPANT  DENSITY  ANP  VENTILATION  RATES 
RESULTING  IN  82°F  DAILY  AVERAGE  EFFECTIVE 
TEMPERATURE  AND  90%  ADEQUACY  FOR  MODEL  I 


maximum 

MAXIMUM 

RANGE 

OF 

NUMBER  OF 

OCCUPANT  DENSITY 

CFM/OCCUPANT  ! 

LOCATION 

OCCUPANTS 

(SQ.FT./OCC) 

B&Hil 

ATLANTA 

28 

21 

229 

PORTLAND 

H 

5 

4 

35 

As  more  extensive  model  testing  is  conducted  for  a  wider  range  of  independent 
variables  and  shelter  characteristics,  more  confidence  can  be  placed  in  the 
results  obtained  from  the  adequacy  assessment  technique. 
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ONIN  AND  9NA*  APE  GIVEN  PE*  ONE  OCCUPANT 

?  **  *  *F  i  g"  re*  *4* !  7*  *  *  *  *SAMPL  E*  *0  F*  DET  A  I*L  ED**ReVu  L  T*S 


SAMPLE  OF  SUMMARY  OF  MAXIMUM  OCCUPANT  DENSITIES 
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Figure  4.9  SAMPLE  OF  SUMMARY  OF  MAXIMUM  AIRFLOW  PER  OCCUPANT 


MAXIMUM  AIK  FLOW  PER  PER  SON  IN  THE  SHELTERS 
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Figure  4.10  SAMPLE  OF  SUMMARY  OF  MINIMUM  AIRFLOW  PER  OCCUPANT 


Reference's  for  Section  4 


4.1  ASHRAE  Applications  Handbook,  Chapter  1C,  "Environmental  Control  For 
Survi val ",  1078. 

4.C  Houohten.  F.  C.,  "Heat  and  Moisture  Losses  From  The  Human  Body  and  Their 
Relation  To  Air  Conditioning  Problems".  ASHVE,  Transactions.  Vol .  35,  1929. 

4.3  ASHRAE  Handbook  of  Fundamentals ,  Chapter  8.  "Physiological  Principles, 
Comfort  and  Health",  1077. 
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Section  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

During  this  experimental  and  analytical  program,  the  following  was  achieved: 

1)  a  specially  designed  low  speed  wind  tunnel  suitable  for  scaled  model 
testing  of  buildings  was  constructed  and  photographic  measurement 
techniques  developed  which  utilized  flow  tracing  of  neutrally  buoyant 
bubbles  through  models 

2)  The  scaled  model  testino  of  three  upgraded  shelter  configurations 
was  completed  establishing  the  ventilation  air  throughput  rate  as  a 
function  of  wind  speed  and  direction. 

3)  An  analytical  flow  model  was  developed  based  upon  the  experimental 
results  and  was  incorporated  into  a  computerized  natural  ventilation 
adequacv  assessment  technique. 

It  is  understood  that  the  air  throughput  rate  is  a  complex  function  of 
many  independent  variables  including  such  dominant  ones  as  wind  speed  and 
direction;  height,  width  and  length  of  the  building;  areas  and  locations  of 
openings  and  internal  resistances.  Certain  of  these  have  been  carefully 
controlled  over  a  range  of  consideration  and  the  resulting  data  is  summarized 
in  Table  6.1.  These  projections  are  made  based  upon  an  assumed  occupant 
density  of  10  sq.  ft.  per  person  and  wind  speeds  in  excess  of  6  mph  (7.3  EPS). 
Extrapolation  of  the  results  to  wind  speeds  lower  than  6  mph  is  not  adviseable 
due  to  the  non-linearity  that  may  exist  in  this  range  and  the  presence  of  other 
dominant  effects  not  yet  considered. 

The  results  summarized  in  Table  6.1  are  very  encouraging  in  the  light  of 
data  reported  hv  f'FMA  (Ref.  6.1)  where  the  per  capita  rate  of  ventilation 
required  to  maintain  an  effective  temperature  of  82 °F  in  an  occupied  space 
with  °(T’  adequacv  durinq  a  normal  year  is  shown  to  range  from  7.6  to  40  CFM 
per  person  depending  upon  climatic  area  in  the  U.S.  All  but  the  southeast 
portion  of  the  U.S.  can  be  satisfied  with  a  rate  as  low  as  20  CFM  per  person. 
Before  these  results  can  be  accepted  with  confidence  however,  the  effects  of 
other  factors,  such  as: 

1)  internal  flow  resistances  including  partition  walls  and  occupants. 

2)  areas  and  locations  of  openings  in  walls, 

l)  wider  selection  of  wind  approach  angles,  and 
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4)  wind  speeds  less  than  5  mph 

should  be  studied  so  that  their  effects  can  be  factored  into  the  overall  results. 
The  existing  test  apparatus  provides  the  capability  to  test  these  parameters  in 
great  enough  detail  and  range  to  allow  accurate  predictions  and  meaningful 
recommendations  to  be  made  with  respect  to  ventilation  of  shelter  buildings. 

The  test  results  also  indicated  an  interesting  effect  that  bermed  openings 
have  on  the  flow  characteristics  of  upgraded  shelters,  i.e.,  larger  shelter 
throughput  rates  are  experienced  when  the  wind  approaches  the  openings  at 
angles  other  than  normal  to  the  opening  (see  results  for  models  I,  II  and  III 
at  45°).  This  effect  should  be  studied  in  greater  depth  and  recommendations 
developed  as  to  how  the  berm  should  be  configured  in  the  neighborhood  of  doors 
and  windows  to  give  best  overall  shelter  throughput  rates. 

In  conclusion,  this  study  had  indicated  that  for  upgraded  shelters, 
ventilation  rates  due  to  wind  forces  are  sizeable,  even  at  low  wind  speeds,  and 
with  the  proper  opening  area  and  location  can  very  probably  be  relied  upon  to 
meet  the  82°F  effective  temperature  and  90°'  adequacy  criteria.  Further 
research  is  necessary  however,  before  the  natural  ventilation  characteristics 
of  upgraded  shelters  are  fully  understood. 
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APPENDIX  B 

VENTILATION  THROUGHPUT  CALCULATION  MODEL 


Appendix  B 

VENTILATION  THROUGHPUT  CALCULATION  MODEL 

Figure  B.l  illustrates  the  application  of  Bernoulli's  equation  to  the 
air  flow  through  a  shelter.  Aj  is  the  flow  area  and  Vj  the  flow  velocity  at 
Section  j  of  the  air  stream  that  eventually  passes  through  the  shelter. 
Applying  Bernoulli's  equation  between  Sections  1  and  2,  we  obtain. 


where,  Pj  and  Pj+j  are  the  static  pressures 
y  the  air  density 

and  H  the  irrecoverable  pressure  head  loss. 


Using  mass  conservation, 

Vi  '  V1V1 


vL  -  A,,,  2 

Therefore,  h^  = 


“\  / 2gh . 

and  Lj  =  Lj+1  =  Aj+lVj+l  =  Aj+1  * 


X 


+U2 


i-HF) 


(Eq.  B.l) 


Denoting  the  smaller  of  the  two  flow  areas  A.  by  A  and  the  larger  area 

J 


by  A  ,  one  gets 


l=a,v 


2gh, 


s  *  A  2  x  s  x 


cv/?gh: 


(Eq.  B.2) 


1  -  (t~) 


Here,  h  is  the  velocity  head  difference  between  the  two  sections,  and 

As  x  At 


x  =  AsV"  ~VST  2  _  A  2 

1  -  (tt)  V\  As 


(Eq.  B.3) 
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Referring  to  Figure  B.l  for  flow  between  sections  1  and  2, 


A2  x  A, 

1  'W-  S 
M1  *2 


and  =  L  =  C^V  2gh^ 


If  h  is  the  velocity  head  difference  between  sections  1  and  7, 

h  _ A  -  *? .  v?  -  A ,  A  -  vf .  A  A .  A  A .  AA .  A:  A 

h  "  2g - Zg - +  “2g - +  ^9  ^9  2g 


i.e.,  h  =  h6  +  h5  +  h4  +  h3  +  h2  +  hj. 


From  continuity. 


clV\‘-  C2'^  h2  "  C3^  h3  =  CW  h4  =  C5"^5  =  C6^  ^6 


C,  2  C.  2  C,  2  C,  2 

h2  =  hl^  *  h3  =  hl^  *  h4  =  hl^  *  h5  =  hl^  ’ 

C1  2 

and  h6  ’  hl(c7> 

b 

rC,2  C,  2  C,  2  C  2  C  2  C,  2  i 
we  get,  h  =  hj  (^-)  +  {^-)  +  (^)  +  (^-)  +  (^-)  +  (^) 


2  2^  1 

In  general ,  h  =  h  C  z 

x  x  i=l  c: 


(Eq.  B.4) 


hx  '  7Z~  x  IfT 
Lx  z  rZ 
i=l  Li 


where  N  =  number  of  flow  resistances,  3  in  Figure  4.4, 

, -  V  2gh 

L  =  Lx  =  CxV  29hx  "  f  2N  1 

vs  y 

i=i  i 
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i.e.,  L  = 


V 


V5qh~ 


2N  1 

1  ?" 
i  =  l  Li 


(Eq.  B.5) 


We  shall  now  show  that  the  specific  area  coefficient  S,  defined  by 
Equation  4.  ,  can  be  expressed  in  terms  of  the  orifice  coefficient  Ci  as, 

J 


S  x  EFe  = 


1 


2N  1 

1  ?" 
i  =  l  S 


,  i  =  1,2 . 2N 


From  4.  , 


1 


-\  /2N  r  "  I  2 
V  I  : 2  -  (e?  I  +  r) 
j=2l - J,L-  z 


,  j  =  2,4, .. . ,2N 


.  -iL 

“j  EFe 


Referring  to  Figure  B.l, 


and 


“2  =  sFe 

a4  zFe 
A6 

a6  “  sFe 


6 

'6 


«2 

(Eq.  B.6) 

a4  EFe 

(Eq.  B.7) 

a6 

(Eq.  B.8) 

EFe  -  Ag  =  EFe(l  -  a^) 

l-a2 

(Eq.  B.9) 
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We  also  had. 


!j'l  =  ApT 


£Fec*2 

B2,L  =  ^  = 


o^EFe 

B2,R  ~  ~K^ 


A4  _  a4EFe 
B4,L  =  AT  ~a0rfe 

3  (J— ) 


U,R  =  A5 


o»L  =  A7 


a^zFe 


agiFe 


a  A  Ire 

(r— 1 

B4,R 


bj  r  » r*-  .  j 

J  ,K  Aj+1 


Al  =  lFe  x  T. 


a2 

A,  =  EFe  x  - - 

J  b2,R 


Bo  n  X  I  80  i  *  Bo  n  X  - 
<:  ,K  »2  h,L  2,K  02 


EFe  X  a. 

Ac  =  - 5 

3  »A  D 


a,  1  -  a? 

B/i  d  x  =  e4,R  (~ J-  ) 


4,R  a. 


6,R  A? 


From  equation  8.3, 


C„.  AI4 

A  A«  -  As 

,2.  AfA2 


A  -  A 
A1  ft2 


(1  -  c^EFe 

s: 


•  A  : 

’  "7 


(1  -  a2)rFe 


(Eq.  B. 10) 


(Eq.  B. 11) 


(Eq.  B.12) 


(Eq.  B. 13) 


(Eq.  B. 14) 


(Eq.  B.15) 
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Substituting  for  Aj  and  A2  from  (B.10)  and  (B.6), 


and 


2 
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Substituting  this  in  Equation  B.5 


~V  2gh 


v: 


(EFe)‘ 


2N 

r 

j-2 


(6j.i  ;  ■*.■>> 

2 

“J 


i .e.  L  = 


V 


2N 

l 

j=2 


2  -  (b! 


EFe,  VTgh 


Bf  d) 


2 

“J 


i.e.  L  = 


S  x  rFe  x  V2gh 


(Eq.  B.22) 


The  equation  derived  above  is  based  on  a  simple  one  dimensional  flow  model. 
In  order  to  account  for  deviations  of  the  real  flow  from  this  simple  model  and 
its  dependence  on  such  variables  as  shelter  orientation,  wind  angle,  shelter 
geometry,  presence  of  earth  berm,  etc.,  the  theoretical  equation  is  modified  by 
incorporating  the  empirical  function  m(v,e,  sc). 

L  -  m  x  S  x  rFe  x  (Eq.  B.23  or  4.2) 
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